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ABSTBACT 

Experimental  data  has  been  obtained  that  predicts  the 
effect  of  roll  angle  upon  the  natural  convection  heat 
transfer  from  an  array  of  cylindrical  spines  or  pin  fins.  A 
correlation  is  presented  that  permits  evaluation  of  the 
Nusselt  Modulus  as  a  function  of  the  Rayleigh  Modulus  and 
the  roll  angle  and  it  shows  that  maximum  heat  transfer  is 
obtained  at  roll  angles  between  20  and  30  degrees. 
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I.  INTRODUCTION 

Electronically  steered  phased  array  antennae  frequently 
exhibit  degraded  performance  due  to  the  lack  of  proper  heat 
transfer  between  the  interior  crystalline  materials  and  the 
environment.  This  causes  higher  than  optimum  temperature 
environments  for  the  crystalline  materials  and  a  subseguent 
loss  of  reliability[  1  ]. 

One   possible  solution   to  such   a   problem  considers   a 
surface  augmentation  through  the  use  of  pin  fins  (or  spines) 
to  assist  in  the  dissipation  of  the  heat. 

The  problem  is  more  severe  when  the  installation  is 
stationary  so  that  the  dissipation  from  the  spines  is  in  the 
natural  convection  mode  and,  because  of  this,  the  focus  of 
this  investigation  is  in  the  area  of  natural  convection  heat 
transfer  and  methods  of  enhancenent  in  this  domain. 

Convective  heat  transfer  from  horizontal  cylinders  has 
been  the  subject  of  extensive  studies  which  were  summarized 
by  Morgan[2j.  In  addition,  studies  of  heat  exchanger  tube 
bundles  are  summarized  by  Incropera  and  and  DeWitt[3]. 
However,  after  a  considerable  search  of  the  literature,  only 
a  small  amount  of  experimental  data  and  analytical  correla- 
tions dealing  with  natural  convection  in  an  array  of  hori- 
zontal cylinders,  has  been  found. 

Intuitively,  one  feels  that  the  orientation  of  the  array 
must  have  an  effect  on  the  heat  transfer  but  apparently  no 
known  experimental  correlations  have  been  formulated.  Thus, 
the  focus  of  this  thesis  is  in  the  effect  of  roll  angle  on 
the  heat  transfer  rates  from  an  array  of  cylindrical  spines. 

Additional  impetus  for  investigations  in  similar  areas 
with  cylindrical  spine  arrays  derives  from  the  electronics 
industry  where  such  cooling  arrangements  are  used  in  the 
external  cooling  of  dense  integrated  circuit  arrangements. 
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Dimension]. ess  correlations  using  dimensionless  parame- 
ters are  usually  developed  for  presentation  of  experimental 
results  and  this  allows  for  possible  application  to  other 
systems.  In  developing  the  solution  for  the  heat  transfer 
across  a  boundary  layer,  it  has  been  shown  by 
Schlichting£4 ],  and  others  [3,5],  that  the  Nusselt  number 
(Nu)  is  a  proper  dimensionless  parameter  for  natural  or 
forced  convection.  In  forced  convection,  the  Nusselt  number 
is  a  function  of  the  Reynolds  number  (Re)  and  the  Prandtl 
number  (Pr)  .  However,  in  pure  natural  convection  where  the 
buoyancy  is  the  only  driving  force,  the  fluid  velocity  is 
determined  entirely  by  the  guantities  contained  in  the 
Grashof  number  (Gr)  which  is  used  instead  of  the  Reynolds 
number  to  establish  dynamic  similarity.  The  Grashof  number 
is  developed  guite  thoroughly  by  Kreith  and  others  [3,5]  and 
quickly  appears  in  a  dimensional  analysis. 

In  addition,  for  gases  the  Prandtl  number  (Pr)  is  nearly 
constant  so  that  the  Nusselt  Modulus  (Nuj)  can  be  considered 
to  be  a  function  of  the  product  of  the  Grashof  number  and 
the  Prandtl  number.  This  product  is  defined  as  the  Rayleigh 
or  number  (Ra)  and  is  discussed  to  further  detail  in  the 
following  paragraphs. 

In  order  to  establish  the  heat  transfer  coefficient  and 
to  develop  a  correlation  of  the  data  with  dimensionless 
parameters,  relationships  for  the  heat  transfer  coefficient 
and  the  Nusselt  number  were  established  and  measurements 
were  taken  which  permitted  its  calculation.  A  summary  of 
the  derivation  is  provided  in  Figure  1. 1  and  a  detailed 
derivation  is  provided  in  Appendix  A.  The  relation  holds 
for  the  condition  that  no  heat  transfer  occurs  at  the  fin 
tip.  Such  a  condition  can  be  effectively  achieved  by  insu- 
lating the  tip. 

After  calculating  the  heat  transfer  coefficient  the 
Nusselt  number  was  then  determined  from  the  relation: 
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where  N u^  is  the  Musselt  number  based  on  the  spine  diameter. 
The  Rayleigh  nunber   vas  determined  from  the   product  of 
the  'irashof  number  and  the   Prandtl  number  which  were  calcu- 
lated from  the  following: 
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Ra  =  GrPr  (1-2) 

vhere  again  the  subscript  "d"  infers  a  dependence  on  the 
spine  diameter  and  where  the  symbols  are  as  defined  in  the 
nomenclature  list. 
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II-  EXPERIMENTAL  APPARATUS 

Figure  2. 1  displays  the  configuration  of  the  experi- 
mental apparatus  and  the  techniques  used  in  data  acquisi- 
tion. The  following  paragraphs  provide  a  more  indepth 
description  of  the  component  parts. 

A.  TEST  DISK 

The  disk,  shown  in  Figures  2-2  and  2.3,  was  fabricated 
entirely  of  2024-T4  aluminum  plate, 6.35  cm  in  thickness.  A 
circular  disk  of  diameter  equal  to  22.86  cm  was  cut  from 
this  plate  and  fitted  with  an  array  of  16  cylindrical  spines 
as  the  heat  transfer  surfaces.  Six  equally  spaced  radial 
holes  were  drilled  around  the  perimeter  of  the  disk  for 
insertion  of  six  200  watt  cartridge  heaters  which  provided 
the  heat  input  necessary  to  elevate  the  temperature  of  the 
spines.  The  heaters  were  powered  by  a  controllable  power 
source  in  order  to  attain  a  steady-state  condition  for  data 
measurements.  The  disk  was  fully  insulated  with  soft  rubber 
refrigeration  insulation  (kx  =  0.15  w/m-K)  to  minimize  heat 
losses  from  the  disk  surface  due  to  both  radiation  and 
convection. 

B.  CYLINDRICAL  SPINES 

The  spines,  shown  in  Figures  2.3  and  2.4  were  similarly 
constructed  of  the  2024-T4  aluninum.  Two  sizes  of  spines 
were  used  in  data  gathering  (0.9525  cm  diameter  x  3.75  cm 
length  and  1.27  cm  diameter  x  5.08  cm  length)  to  provide  a 
wider  range  of  Rayleigh  Numbers.  Extension  pieces,  also 
made  of  the  same  aluminum,  were  necessary  to  extend  the 
actual   base      of    the    heat      transfer   surface       (the    cylindrical 
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Figure    2.1         Sketch    of    the   Test    Apparatus    and   Configuration 
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n  I, inimical 

SPIN!'. 


Figure  2.2    Test  Disk  and  Spine  Assembly 

spine)    to   the   outer  level   (exterior   surface)    of   the 

insulation . 

C.   THERMOCOUPLES 

All  thermocouples  were  fabricated  from  copper-constantin 
(Type-T)  of  AWG  No.  30  wire.  They  were  butt  welded  with  a 
spherical  bead  of  approximately  0.05  cm  diameter  and 
smaller.  All  interconnecting  cables  were  also  of  copper- 
constantin  in  order  to  minimize  the  uncertainty  and  to 
preclude,  as  much  as  possible,  spurious  thermoelectric 
effects.  All  thermocouples  were  connected  to  the  ACUREX 
Autodata  Nine  data  acquisition  system  for  temperature 
recording. 
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Figure    2.3         Sketch   of    Front    View    of   Test    Disk 

D.       ACOREX    ADTODATA    NINE 

The  thermocouple  modules  in  the  data  acquisition 
system's  thermocouple  modules  include  an  ice-point  compensa- 
tion network  and  provide  automatic  data  scanning  and 
printout      for    up      to    100      channels      of    input.  Due    to      the 

internal  ice-point  compensation  networks,  no  external  ice 
baths  or  reference  junctions  were  required.  In-place  cali- 
bration checks  were  made  at  various  temperatures  over  the 
expect€d  temperature  range  in  crder  to  develop  temperature 
correction  factors  for  each  thermocouple,  its  connecting 
cabling  and  the  Autodata  Nine  circuitry  as  a  unit.  in  addi- 
tion, the  complete  disk  and  spine  array  was  immersed  in  an 
ice-bath   in    order   to   establish    cne    of    the    check    points. 
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figure    2.4         lest    Disk    and   Spine    Arrangement 

E.       ENVIRONMENTAL    CHAMBER 

The  test  system  was  housed  in  an  environmental  chamber 
in  which  the  ambient  temperatures  could  be  maintained  over  a 
range   between    35    and    135    degrees   Fahrenheit. 


18 


III.  EXPERIMENTAL  PROCEDURE 

A.   DATA  GATHERING 

The  experimental  procedure  was  to  set  the  disk  at  the 
desired  roll  angle  using  a  liguid  bubble  level  indicator  and 
protractor  to  measure  the  angle.  The  temperature  control 
thermostat  of  the  environmental  chamber  was  then  set  to  the 
desired  temperature.  Steady  state  environmental  conditions 
were  attained  in  approximately  six  to  eight  hours.  Next  the 
heaters  imbedded  in  the  disk  were  energized  to  elevate  the 
temperature  of  the  disk  and  base  of  the  spine  to  a  tempera- 
ture chosen  for  a  predetermined  Rayleigh  number.  The  input 
voltage  for  the  heaters  was  varied  and  the  temperatures 
monitored  until  a  steady  state  condition  was  reached.  This 
was  indicated  by  constant  input  power  and  constant  tempera- 
tures. The  time  to  attain  a  total  steady-state  condition 
was  approximately  three  to  five  hours  in  addition  to  the 
time  required  for  the  chamber  tc  stabilize. 

A  steady  state  condition  was  assumed  to  exist  when  no 
change  in  temperature  occurred  for  four  consecutive  tempera- 
ture readings  taken  at  fifteen  minute  intervals.  These 
temperatures  were  then  used  for  the  input  data  for  the 
calculations.   Typical  dataruns  are  presented  in  Appendix  B. 

After  completion  of  data  taking  at  a  set  condition, 
changes  were  made  in  the  roll  angle,  the  temperature  of  the 
system  or  temperature  of  the  environment,  and  the  process  of 
attaining  and  maintaining  a  new  condition  of  steady  state 
was  again  started. 
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B.   DATA  SEDUCTION 

Temperature  data  printed  by  the  Autodata  Nine  during  the 
individual  runs  was  used  as  input  parameters  (corrected  from 
the  calibration  curve  for  each  thermocouple)  for  calculation 
of  the  desired  parameters  for  plotting  and  evaluation. 

All  desired  output  parameters  were  calculated  by  the 
computer  using  appropriate  relationships  for  the  thermophy- 
sical  properties  [3]  as  functions  of  temperature  and  meas- 
ured temperatures  as  input  paraueters.  The  computer  program 
used  for  the  calculations  is  provided  in  Appendix  E. 

Graphical  display  of  the  data  was  also  utilized  to  aid 
in  data  analyses.  These  curves  are  shown  as  Figures  B.  1 
through  B.7. 
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IV.     RESULTS    AND   CORE ELATION 

A.  RESULTS 

During  each  data  collection  run  the  temperatures  at  the 
base  and  tip  of  each  spine,  as  well  as  four  environmental 
temperatures  were  recorded.  These  temperatures  were  used  as 
inputs  for  calculation  of  the  various  output  parameters  for 
each  spine.  These  sixteen  calculations  were  then  averaged 
to  provide  the  data  point  for  that  set  of  conditions.  Forty 
seven  data  points  were  obtained.  Typical  data  runs  are 
included  as  Appendix  C. 

At  a  specified  angle,  the  natural  logarithm  of  the 
Nusselt  number  was  plotted  as  a  function  of  the  natural 
logarithm  of  the  Rayleigh  number  for  each  of  the  data  points 
at  a  specified  angle.  These  jlots  are  included  as  Figures 
B. 1  through  B.5  in  Appendix  B.  This  graphical  display  was 
then  used  in  order  to  facilitate  correlation  of  the  test 
data . 

The  plotted  results  indicated  that  a  small  dependence  on 
roll  angle  did  exist.  Figures  B.6  and  B.7  were  plotted  to 
show  these  effects  more  specifically.  Figure  B.7  was 
plotted  only  for  Rayleigh  nunbers  between  579  and  586 
because  of  the  narrowness  of  the  overall  Rayleigh  number 
range  obtained. 

B.  CORRELATION 

For  each  of  the  curves  showing  log  (Nu)  as  a  function  of 
log  (Ra) ,  a  least  squares  curve  fit  routine  was  used  and  a 
linear  relation  established  to  approximate  these  curves 
following  a  method  presented  by  Kraus  and  Bar-Cohen[ 1  ].  A 
correlating  equation  was  obtained  for  each  of  the  roll 
angles  (0),  where  Z  [$)    defines  a  function  of  the  roll  angle. 
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The  values  or  Z  (j^)  are  tabulated  ia  Table    1 


TABLE  1 
Enhancement  Factor 


Roll 

A  i 

igle(0) 

0° 

10° 

20° 

30° 

45° 

Enhancement 

Factor  (Z(0)) 

0 

.5916 

0 

.5986 

0 

.6075 

0 

.5990 

0 

.5928 

The  values  of  Nusselt  numbers  obtained  from  the  experi- 
mental results  are  a  function  of  a  giver.  Bdyleigh  number  and 
roll  angle  and  may  re  correlated  via: 


NUj  =  Z(0)Ra 
d 


0.233 


(U.1) 


However,  observe  that  the  dependence  on  the  roll  angle  is 
weak  (Table  1)  and  no  mathematical  relationship  was  deter- 
mined for  the  values  or  Z {$) . 

Results  were  also  compared  with  several  previously 
published  correlations  for  the  single  horizontal  cylinder. 
Churchill  and  Chu£6]  suggested  an  equation  ot  the  form 


Nu 


0.60  + 


0 .387  Ra 


(1/6) 


1  + 


0.559 
Pr 


(9/16) 


(8/27) 


22 


which  correlated  experimental  results  of  several  researchers 
very  well. 

The  values  obtained  in  this  experimental  program  agreed 
favorably  with  these  predictions  with  an  average  of  19% 
deviation  throughtout  the  range  investigated.  Noting  that 
the  effect  of  heat  transfer  by  radiation  were  also  present, 
the  deviation  is  attributed  to  radiation  losses  primarily. 
These  losses  would  appear  to  net  have  a  detrimental  effect 
on  the  observed  roll  angle  dependence. 
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V.  CONCLUSION  AND  EECOMMENDATIONS 

A.  CONCLUSION 

The  results  of  the  experimental  program  indicated  that  a 
dependence  on  roll  angle  did  exist.  This  dependence  is 
relatively  minor  (i.e.,  on  the  order  of  3%),  however  for 
temperatures  in  the  range  of  200  degrees  Fahrenheit  this 
could  mean  a  difference  of  5  to  6  degrees  in  the  surface 
temperature  of  an  operating  system.  Therefore,  further 
investigation  at  higher  Rayleigh  numbers  may  be  deemed 
appropriate. 

B.  RECOMMENDATIONS 

The  following  recommendations  are  provided  for  possible 
follow  on  projects  of  a  similar  nature. 

1.  Conduct  the  experimental  prcgram  with  the  goal  of  deter- 
mining and  removing  the  radiation  heat  transfer. 

2.  Investigate  other  devices  and  methods  of  temperature 
measurement  in  order  to  remove  more  of  the  uncertainty 
associated  with  the  use  of  thermocouples. 

3.  Conduct  a  similar  experiment  for  the  pitch  angle. 

4.  Attempt  a  similar  program  with  water  as  the  convecting 

medium. 
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APPEND 21  4 
DEFINITION  OP  THE  HEAT  TRANSFER  COEFFICIENT 

Fins  and  spines  (pin  fins)  of  various  geometries  and 
thermal  conductivities  respond  differently  to  identical  and 
uniform  heat  sources  and  sinks.  Similarly,  there  are 
numerous  ways  in  which  the  temperatures  and  heat  transfer 
coefficients  of  sources  and  sicks  may  vary.  Important  to 
the  analysis  of  fin  geometries  are  the  constraints  or 
assumptions  which  are  employed  to  define  and  limit  the 
problem  and  often  to  simplify  its  solution.  The  analysis  of 
fins  and  spines  employ  the  assumptions  of  Murray[7]  and 
Gardner[8].   These  limiting  assumptions  are: 

1.  The  heat  flow  in  the  fin  and  its  temperatures  remain 
constant  with  time. 

2.  The  fin  material  is  homogeneous,  its  thermal  conduc- 
tivity is  the  same  in  all  directions  and  it  remains 
constant. 

3.  The  heat  transfer  coefficient  between  the  fin  and 
the  surrounding  medium  is  constant  and  uniform  over 
the  entire  surface  of  the  fin  (spine) . 

4.  The  temperature  of  the  medium  surrounding  the  fin  is 
uniform. 

5.  The  fin  thickness  is  so  snail  compared  with  its  height 
that  temperature  gradients  across  the  fin  thickness 
may  be  neglected. 

6.  The  temperature  at  the  base  of  the  fin  is  uniform. 

7.  There  is  no  contact  resistance  where  the  base  of  the 
fin  joins  the  prime  surface. 
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8.  There  are  no  heat  sources  within  the  fin  itself. 

9.  The  heat  transferred  through  the  outermost  edge  of 
the  fin  is  negligible  compared  with  that  leaving  the 
fin  through  its  lateral  surface. 

10.  Heat  transfer  to  or  from  the  fin  is  proportional  to 
the  temperature  excess  between  the  fin  and  the  sur- 
rounding medium. 

In  addition  to  proposing  a  profile  function  for  longitu- 
dinal fins  Gardner[8]  proposed  a  profile  function  for 
spines: 


(-2. 


f2(x)  = 


^(-r-r 


(A.1) 


With  the  appropriate  value  of  n,  eguation  A.1  may  be  used 
for  the  development  of  a  generalized  differential  equation 
for  spines. 

Figure  A.1  shows  a  spine  of  arbitrary  profile.  It  can 
be  seen  that  the  spine  cross  section  normal  to  the  flow  of 
heat,  the  confining  profile,  and  the  perimeter  of  the  spine 
are  perfectly  arbitrary  functiors  of  the  distance  x  from  the 
tip  of  the  spine.  A  differential  equation  for  the  tempera- 
ture excess,  0  =  t  -  t  ,  may  re  written  by  considering  the 
heat  conduction  into  and  out  of  the  element  dx  through  the 
cross  section  f  (x)  : 


dq   =  k 


dx 


fl(x)  3 


I]dx 
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x  =  i> 


Jt-0 


Figure  A. 1    Spine  of  Arbitrary  Profile 

The  difference  in  heat  flow  into  and  out  of  the  element 
must  equal  that  dissipated  by  the  surface  of  the  spine.  If 
dissipation  occurs  by  convection  and  h  is  the  convection 
coefficient 

dq  =  h  fgCx)  0  dx 


where  f3(x)  defines  a  perimeter  function  P  (:<)  which  depends 
upon  the  distance  x  from  the  origin  of  the  coordinate 
system.   When  convection  and  conduction  are  equated, 


kd 


x  I  2    dx 


hf 3(x)9 


Upon   rearrangement   the  generalized   differential   equation 
becomes 


f  (x) 


d2Q 


dx' 


df  (x)  dO 


dx     dx     k 


-f3(x)  =  0 


(A. 2) 
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T h e  relationship  between  f ,  ( x )  and  f  ( x )     i 


f^x)  =  7T  [f2(x)l 


2 


where  f2(K)  w-i-3  defined  by  equation  A.I.    Thus  equation  A.  2 
may  be  written 


[f2(x)]2  d20    d  rf.(x)''2 


— 2"  +  — _ 
dx     dx 


dQ     2h  f2(x)9         (A.3) 
dx     k 


where 


f3(x)  =  2rrf2(x) 


Equation  A. 3  is  a  second-order  differential  equation 
with  variable  coefficients  except  when  the  spine  cross 
section  normal  to  heat  flew  is  constant.  It  may  be  solved 
by  termwise  comparison  with  the  generalized  Hossel  equation. 

For  the  cylindrical  spine,  shown  in  Figure  A. 2 


V*)  =  -f 


and 


a«[f2<*>] 


=  0 


For  this   case  the  exponent  in   equation  A.1  has  a   value  of 
zero  correspond  in j  to  n  =   V2.    Actually  £    in  equation  A.1 
is  replaced  by  d,  the  spine  diameter.   Using  this  value  of  f2 
the  general   differential  equation   for  the   temperature 
excess  becomes 
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Figure   A. 2         Cylindrical    Spine 


d29 


dx 


-    m   e    =    o 


(A. 4) 


where    m    =     (4h/kdJ*  . 

The    general    solution    to    equation    A. 4    is 


9(x)    =    C,emx    +      C    e"mx 


(A. 5) 


where   the   arbitrary  constants   are   determined   from   the 
boundary  conditions 


9Cx=b) 
9(x=0) 
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=  9. 


(A. 6a) 

(A. 6b) 


This  then  leads  directly  to    the  equation  for  the  temperature 
excess  profile 


(x)  =  9 


cosh  mx 
b  cosh  mb 


(A. 7) 


and  the  heat  flow  entering  the  fin  (spine)  base  is 


qb  =  -J— kd  m9   tanh(mb) 


(A. 3) 


The  value  of  h  may  be  determined  from  equation  (A. 7)  by 
taking  temperature  readings  at  the  tip  and  the  base  which 
are  precisely  "b"  meters  apart. 

Thus  with  base  and  tip  temperature  excesses, 

g  (x=b)  =  ebe-(x=aj  =  e^ 

It  is  easy  to  see  from  equation  (A. 7)  that 


9, 


9  = 


a"  cosh  mb 


or  that 


9, 


cosh(mb)  = 


from  which 


mb  =  arccosh(9  /9  ) 
d   a 


=  a 


(A.  9) 


Then  with   m  =  (4h/kdj2.   one  may  then  solve  for  h 


h  =  — (b} 


30 


APPENDIX  B 
GRAPHICAL  DISPLAY  OF  DATA 
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APPENDIX  C 
TYPICAL  DATA  RUNS 
TABLE  C.l 
Typical  Results  (579  <  Ra  <  586) 


T  = 
s 

296  .  89 

0  = 

:  0° 

h 

Pin 

TbC< 

312. 

'K) 

,48 

T  (°K) 
o 

Ra 

Nu 

Nu  . 
ch 

2.3333 

(w/r 

6  , 

2 

1 

312, 

.25 

579  , 

,40 

2 

.2339 

.162 

2 

312  , 

,  32 

312, 

.05 

579  , 

,  31 

2, 

.5850 

2  . 

.  3283 

7  , 

.128 

3 

312  , 

,50 

312 

.23 

579  , 

,40 

2 

.9536 

2  , 

.3333 

7  , 

,043 

it 

312  , 

.  32 

312 

.03 

572  , 

.94 

2, 

.  8039 

2  , 

.  3280 

7  , 

.731 

t 

o 

312, 

,  37 

312 

.08 

574  , 

.69 

2 

.7417 

2  , 

.  3295 

7  , 

.560 

h 

312  . 

.59 

312 

.  33 

582  , 

,43 

2 

.4842 

2  , 

.  3358 

6 

.853 

7 

312  , 

,61 

312 

.  3  5 

583, 

.17 

2 

.4806 

2  , 

.  3364 

6 

.843 

b 

312  , 

,64 

312 

.  3  4 

583, 

,  84 

2 

.  7984 

2  , 

.  3367 

7 

.720 

9 

312  , 

.HI 

312 

.16 

576  , 

,73 

2 

.4066 

2  , 

.3311 

6, 

.637 

10 

312  , 

,  14  M 

312 

.16 

577, 

,28 

2 

.  7274 

2  , 

.3316 

7, 

.522 

]  I 

312  , 

,58 

312 

.  31 

581, 

,97 

2 

.6458 

2  , 

.  3354 

7  , 

.298 

]  2 

312  , 

,66 

312 

.  38 

584  , 

,27 

2 

.  5824 

2  , 

.  3373 

7, 

.124 

]  3 

312. 

.60 

312 

.  32 

582  , 

.43 

2 

.6995 

2  , 

,  3358 

7  , 

.446 

14 

312  , 

.  37 

312 

.  14 

5  75  , 

.71 

2 

.2505 

2  , 

,3303 

6  , 

.206 

]  5 

3  1  2  . 

,  5  4 

312 

.25 

5  8  0, 

.  31 

2 

.  7625 

2  , 

.  3340 

7  , 

.620 

If. 

312, 

,48 

312 

.19 

578, 

.48 

2 

.7727 

2 

.3326 

7  , 

.647 

T  =  2  9  6.93   0=lu° 

1  312.76  312.47  586.28  2.7304  2.3389  7.534 

2  312.5  9  312.33  5  81.14  2.5429  2.33  4  7  7.015 

3  312.78  312.51  587.20  2.5123  2.3396  6.932 

4  312.59  312.31  580.77  2.7581  2.3344  7.609 

5  312.6  4  312.36  582.51  2.6972  2.3358  7.441 

6  312.81  312.55  588.37  2.4534  2.3406  6.770 

7  312.89  312.63  590.95  2.4434  2.3426  6.743 
H  312.86  312.57  589.49  2.7638  2.3415  7.627 
9  312.63  312.38  582.70  2.3763  2.3360  6.556 

10  312.67  312.3  8  583.25  2.6932  2.3364  7.430 

11  312.86  312.58  589.76  2.6057  2.3417  7.191 

12  312.8  7  312.60  588.37  2.5903  2.3420  7.038 

13  312.82  312.54  590.21  2.6658  2.3406  7.356 

14  312.59  312.31  580.77  2.7581  2.3344  7.609 
16  312.7  6  312.47  5  80.25  2.7304  2.3389  7.534 
16  312.71  312.4  2  584.44  2.7380  2.3374  7.554 
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TABLE  C.l 
Typical  Results  (579  <  Ra  <  586)  (cont'd.) 


T   =296.93 


7 

H 

9 

10 

]  1 
12 
1  3 
1  'I 

15 
16 


312 
312 
312 
312 
312 
312 
312 
312 
312 
312 
312 
312 
31  2 
312 
312 
312 


71 
54 
72 
5  4 
59 
7  6 

7  8 
81 
58 
56 
75 

8  2 
77 
5  4 
59 
5  9 


=  20 

312 
312 
312 
312 
312 
312 
312 
312 
312 
312 
312 
312 
312 
312 
312 
312 


42 
27 
46 
25 
31 

4  4 
52 
51 
33 
27 
47 

5  4 
48 
2  5 
31 
31 


584 
579 
585 
578 
580 
565 
587 
587 
580 
579 
586 
588 
586 
578 
5  80 
580 


44 
29 
36 
92 
68 
84 
28 
65 
86 
57 
09 
3  7 
56 
92 
77 
77 


7380 
5523 
5216 
7684 
7047 
9965 
4611 
7738 
3851 
7130 
6246 
5597 
6733 
7680 
7584 
7584 


3374 
3332 
3381 
3329 
3343 
3384 
3397 
3400 
3345 
3343 
3387 
3406 
3391 
3329 
3344 
3  34  4 


554 
040 
958 
636 
461 
268 
791 
654 
579 
484 
242 
063 
376 
636 
609 
609 


T 


=  296 .94 


0  =  30° 


7 

8 

9 

10 

11 

12 

1  3 
14 
]  5 
lb 


312 
312 
31  2 
312 
312 
312 
312 
312 
312 
312 
312 
312 
312 
312 
312 
312 


65 
54 
67 
5  4 
5  3 
76 
7  8 


70 
69 
5  2 
7  7 
48 
59 


312 
312 
312 
312 
312 
312 
312 
312 
312 
31  2 
3J2 
312 
312 
312 
312 
312 


36 
27 
40 
25 
31 
44 
52 
27 
27 
51 
42 
24 
4  8 
14 
31 
31 


582 
579 
583 
578 
579 
58  5 
587 
579 
579 
586 
584 
588 
586 
575 
5  80 
580 


81 
29 
52 
92 
7  6 
64 
28 
94 
97 
73 
27 
37 
56 
25 
77 
77 


7480 

5523 
5332 
7684 
1768 
99  6  5 
4  611 
9250 
7130 
2509 
6842 
6599 
6733 
7890 
7581 
7581 


3359 
3332 
3366 
3329 
3336 
3384 
3397 
3337 
3334 
3392 
3373 
3406 
3391 
3299 
3344 
3344 


7.581 
7  .040 
6  .989 
7.636 


8.268 

6.791 
8  .068 
7  .484 
6.211 
7.268 
7  .063 
7.  376 
7.692 
7.609 
7  .609 
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TABLE  C.l 
Typical  Results  (579  <  Ra  <  586)  (cont'd.) 


T   =  296.96   0 


45' 


1 

312, 

.65 

312 

.42 

582, 

.12 

2, 

.21640 

2, 

.  33550 

6  , 

.115 

2 

312  , 

.48 

312 

.22 

576, 

.06 

2, 

.56722 

2, 

.  33057 

7 

.081 

3 

312  , 

.67 

312 

.40 

582, 

.12 

2 

.53847 

2  , 

.33550 

7  , 

.004 

4 

312, 

,M8 

312, 

.19 

575, 

.69 

2 

.78452 

2, 

.  33024 

7  , 

.681 

5 

312, 

.  53 
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.25 

577  , 

.44 
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.72293 

2, 

,  33169 

7, 

.511 

6 

312, 

.  7  6 

312 

.49 

585, 

.15 

2 

.46985 

2  , 

.33796 

6 

.  815 

7 

312, 

.  78  ' 

312, 

.  52 

5  85, 

,89 

2 

.46626 

2  , 

.33856 

6  , 

.806 

8 

312, 

.  7  5 

312, 

.51 

585  , 

.  35 

2 

.25553 

2  , 

.33812 

6  , 

.2  24 

9 

3  1 2 

.58 

312 

.  33 

5  79  , 

.47 

2 

.39009 

2, 

, 33335 

6 

.594 

10 

312 

.50 

3 1 2 

.22 

576  , 

.  3  4 

2 

.72653 

2  , 

.33079 

7, 

.521 

1  1 

312, 

.69 

312, 

.42 

582  , 

.86 

2, 

.63975 

2  , 

.  33611 

7, 

.283 

1  '/ 

312, 

.  82 

312 

.54 

5  8  6, 

.99 

2 

.56502 

2  , 

.33945 

7 

.078 

J  3 

3.1  2 

.71 

312 

.43 

583, 

.  31 

2 

.69100 

n 

.  336M7 

7 

.425 

]4 

312  , 

.43 

312 

.14 

573, 

.  85 

2, 

.79496 

2  , 

,32875 

7  , 

.709 

IS 

312  , 

.59 

312  , 

.  31 

579  , 

.  38 

2 

.76402 

2  , 

.  33328 

7 

.625 

16 

312 

.59 

312 

.  31 

579  , 

.38 

2. 

.76402 

2  , 

.  33328 

7 

.625 
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TABLE  C. 2 


Typical  Results  (1305  <  Ra  <  1380) 


T  =317 
s 

.83 

0  = 

0° 

h 

Pin 

Tb(< 
338 

49 

T  (°K) 
o 

338.04 

Ra 

Nu 

Nu  . 
ch 

2.76085 

(w/ 

6 

m 

1 

1304 

18 

2 

.97416 

553 

2 

338 

49 

337 

.99 

1302 

95 

3 

.27379 

2 

76029 

7 

213 

3 

338 

.52 

338 

.06 

1305 

56 

3 

07919 

2 

76146 

6 

784 

4 

338 

.44 

337 

.04 

1300 

20 

3 

28048 

2 

75905 

7 

227 

5 

338 

.53 

338 

.05 

1305 

56 

3 

.15390 

2 

.76146 

6 

949 

6 

338 

.  u  G 

338 

.02 

1303 

.57 

3 

.12388 

2 

.76057 

6 

.882 

7 

338 

.52 

338 

.05 

1305 

40 

3 

.11699 

2 

.76139 

6 

.868 

8 

338 

.50 

338 

.01 

1303 

72 

3 

.23238 

2 

.76064 

7 

.122 

9 

338 

.53 

338 

.08 

1306 

47 

3 

.00352 

2 

.76188 

6 

.618 

10 

338 

.53 

338 

.08 

130  6 

31 

2 

.96674 

2 

76180 

6 

537 

11 

338 

6  7 

338 

.21 

1313 

81 

2 

.98347 

0 

.76517 

6 

575 

12 

338 

6  3 

338 

17 

1311 

6  6 

3 

.06200 

2 

.76421 

6 

74  7 

13 

338 

50 

333 

.01 

1303 

73 

3 

.23233 

2 

.76064 

7 

122 

14 

3  38 

51 

338 

0  4 

1304 

.78 

3 

04564 

2 

76111 

6 

.710 

1  5 

33  3 

.55 

338 

07 

1306 

63 

3 

.18820 

2 

.76195 

7 

025 

1.6 

338 

51 

3  38 

01 

13  0  4 

02 

3 

3  0  6  4  2 

2 

76077 

7 

285 

T  =314 
s 

.72 

0  = 

10° 

1 

334 

.  82 

334 

.  35 

1326 

67 

3 

.23972 

2 

.77102 

7 

.075 

2 

3  3U 

.  87 

3  34 

.42 

13  30 

.20 

3 

.07632 

2 

.77259 

6 

.718 

3 

334 

.81 

334 

.33 

1325 

.  54 

3 

.28189 

2 

.77052 

7 

.166 

4 

334 

.  79 

334 

.  35 

1325 

88 

3 

.04951 

o 

.77067 

6 

.659 

5 

3  34 

.76 

334 

.32 

1323 

93 

3 

.97872 

2 

.  76980 

6 

.504 

6 

3  3  4 

.82 

334 

.  33 

1327 

.63 

3 

.00764 

2 

.77145 

6 

.568 

7 

3  34 

.  75 

3  34 

.28 

1322 

.49 

3 

.24195 

2 

.76916 

7 

.101 

8 

3  34 

.88 

334 

.39 

1329 

.72 

3 

.  34543 

2 

.77238 

7 

.306 

9 

3  34 

.87 

3  34 

.  4  2 

1330 

.06 

3 

.03809 

2 

.77253 

6 

.635 

10 

334 

.  81 

334 

.  37 

1326 

.  84 

3 

.05689 

2 

.77110 

6 

.654 

11 

3  34 

.  8  8 

334 

.41 

1329 

.89 

3 

.23042 

2 

.77245 

7 

.055 

12 

3  34 

.92 

334 

.41 

1331 

.16 

3 

.45717 

2 

.77302 

7 

.  550 

1  3 

3  34 

.95 

3  34 

.49 

1334 

.06 

3 

.06572 

2 

. 774  30 

6 

.696 

14 

334 

.  80 

334 

.  3  4 

1325 

.88 

3 

.12721 

2 

.77067 

6 

.829 

1!, 

3  34 

.93 

334 

.47 

13  33 

.  24 

3 

.18207 

2 

. 77394 

6 

.950 

16 

3  34 

.  80 

3  34 

.33 

13  2  5 

.54 

3 

.20581 

2 

. 77052 

7 

.000 

K) 
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TABLE  C.2 
Typical  Results  (1305  <  Ra  <  1380)  (Cont.) 


T  =311.24   0 


10° 


1 
2 

3 
4 

r 

6 

7 

8 

9 

10 

11 

]2 

13 

14 

15 

1C 


330 

3  30 
3  30 
330 
3  30 
3  30 
3  30 
3  30 
3  30 
3  30 
3  3d 
330 
3  30 
3  30 
330 
3  30 


59 
'i  8 
4  8 

4  3 

5  9 
5  4 
4  3 
59 
54 

4  3 

5  9 

6  0 
53 
4  8 
59 
37 


330 
330 
330 
3  2  9 
330 
3  30 
330 
3  30 
330 
3  29 
3  30 
330 
330 
330 
330 
329 


16 

07 
04 
98 
14 
07 
00 
16 
09 
98 
14 
14 
08 
05 
16 
92 


1347 
1340 
1339 
13  36 
1346 
1342 
1337 
1347 
1343 
]  336 
1346 
134  6 
1342 
1340 
1346 
1333 


11 
81 
97 
56 
6  0 
69 
07 
11 
38 
5  6 
60 
60 
87 
30 
93 
00 


11643 
01289 
21805 
22605 
23823 
37037 
10318 
116  5  3 
20669 
22605 
23823 
31969 
25068 
13528 
15756 
27927 


78019 
77741 
77704 
77553 
77997 
77824 
77576 
78019 
77855 
77553 
77997 
7  7  9  9  7 
77832 
77719 
78012 
77396 


736 
512 
955 
972 
000 
285 
706 
736 
931 
972 
000 
176 
026 
776 
825 
086 


T  =310 
s 


65   0  =  20' 


1 

8 

9 

10 

11 
12 
1  3 
14 
15 
16 


330 
3  30 
3  30 
3  30 
3  3  (J 
330 
3  30 
3  3  f) 
3  30 
3  30 
330 
330 
330 
330 
330 
330 


32 
32 
26 
2  6 
32 
26 

2  6 

3  7 

4  3 
32 
4  3 
43 
42 
32 
4  3 
26 


329 

329 
329 
329 
329 
329 
3  29 
3  29 
3  29 
329 
329 
329 
3*2  9 
329 
329 
329 


89 
82 
82 
87 
80 
78 
88 
9  3 
87 
98 
98 
93 
82 
98 
82 


1377 
1377 
1373 
1373 
1376 
1372 
1372 
1378 
1381 
13  76 
1383 
1383 
1381 
1375 
1383 
1373 


21 
54 
47 
47 
70 
95 
26 
91 
96 
88 
67 
85 
82 
17 
67 
4  7 


06971 
98876 
16018 
16018 
18968 
28058 
44194 
46225 
53134 
14908 
13060 
16992 
49360 
55229 
13060 
16018 


79335 
79350 
79173 
79173 
79313 

79151 
79121 
79409 
79541 
79321 
79614 
69622 
79535 
79247 
79614 
79173 


628 
453 
823 
823 
887 
082 
431 
476 
625 
799 
76  0 
845 
544 
669 
760 
823 
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TABLE  C.2 
Typical  Results  (1305  <  Ra  <  1380)  CCont.) 


T  =303.00   0 


30° 


] 

319 

.  53 

319 

.14 

1317 

29 

3 

31536 

2 

76727 

6 

983 

2 

319 

.  53 

319 

.18 

1318 

68 

3 

07108 

2 

76789 

6 

468 

3 

319 

.56 

319 

.18 

1319 

86 

3 

2  5  8  9  4 

2. 

76841 

6. 

864 

'i 

319 

.48 

319 

.09 

1313 

92 

3 

37431 

2 

76576 

7 

106 

c 
J 

319 

.55 

319 

.18 

1319 

48 

3 

16447 

2 

76824 

6 

665 

6 

319 

.47 

319 

.13 

1314 

72 

2 

98398 

2 

76611 

6 

285 

7 

319 

.51 

319 

.12 

1315 

73 

3 

.31993 

2 

76657 

6 

992 

8 

319 

.  52 

319 

.14 

1316 

70 

3 

.26802 

2 

76700 

6 

883 

9 

319 

.  55 

319 

.17 

1319 

06 

3 

26120 

0 

L 

76806 

6 

869 

10 

319 

.  4  9 

319 

.12 

1315 

10 

3 

27261 

2 

76629 

6 

892 

11 

319 

.52 

319 

.17 

1317 

88 

n 

.97346 

9 

L 

7  6  7  5  3 

6 

263 

12 

319 

.  44 

319 

.09 

1312 

56 

3 

.03882 

2 

76515 

6 

400 

13 

319 

.  56 

319 

.22 

1321 

.04 

2 

96517 

2 

.76894 

6 

246 

14 

319 

.  54 

391 

.17 

1318 

47 

3 

.21405 

2 

.76779 

6 

.770 

i£ 

319 

.47 

319 

.08 

1312 

.95 

3 

,42424 

2 

76532 

7 

211 

16 

319 

.  57 

319 

.22 

1321 

46 

3 

.06151 

2 

76913 

6 

.449 

T  = 

s 

:299  .02 

0  = 

30° 

1 

31 '( 

.53 

314 

.19 

1320 

.50 

3 

.10  591 

2 

.76886 

6 

.462 

2 

3  1  'i 

.53 

314 

.20 

1320 

.91 

3 

.10715 

2 

.76904 

6 

.4  64 

3 

3  1  4 

.  5  7 

314 

.23 

1  323 

.44 

3 

.  20306 

2 

.77017 

6 

.664 

M 

314 

.  57 

314 

.24 

1323 

.66 

3 

.04478 

2 

.77027 

6 

.  335 

i 

31 'J 

.  5  6 

314 

.23 

1322 

81 

3 

.04922 

2 

.76989 

6 

.  344 

f, 

314 

.  4  7 

314 

.  14 

1316 

.47 

3 

.06411 

2 

.  76706 

6 

.  374 

7 

314 

.  5  2 

314 

.18 

1319 

.64 

3 

.10825 

2 

.76848 

6 

.466 

H 

314 

.  5  6 

314 

.23 

1323 

.03 

3 

.09906 

2 

.  76999 

6 

.  448 

'J 

3]  4 

.61 

314 

.29 

1327 

.46 

2 

.93205 

2 

.77196 

6 

.101 

10 

314 

.  50 

314 

.18 

1319 

.00 

2 

.9540s 

2 

.76819 

6 

.146 

1  1 

314 

.  5  6 

314 

.23 

1323 

.25 

3 

.04584 

2 

.  77009 

6 

.  337 

12 

314 

.51 

314 

.14 

1317 

.  51 

3 

.42918 

2 

.76752 

7 

.134 

1  3 

314 

.49 

314 

.12 

1316 

.02 

3 

.48656 

2 

.76686 

7 

.253 

14 

314 

.41 

314 

.15 

1318 

.15 

3 

.  37434 

2 

.76781 

7 

.020 

15 

3]  4 

.47 

314 

.10 

1314 

.57 

3 

4  3  8  0  9 

2 

7  6  6  21 

7 

152 

1G 

314 

.  53 

314 

.20 

1320 

.91 

3 

.10715 

2 

.76904 

6 

.464 

44 


TABLE  C.2 
Typical  Results  (1305  <  Ra  <  1380)  ( Cont . ) 


T    =301.86 
s 


1 
2 
3 

4 

r 

.j 

6 
7 

8 
9 
10 
11 
12 
13 
14 
15 
lb 


318 
318 
31  8 
318 
319 
31  8 
319 
318 
31  8 
31  8 
318 
319 
318 
318 
318 
318 


91 
87 
9G 
9  6 

00 
9  7 
02 
95 
89 
9  4 
95 
01 
99 
96 
91 
98 


30  .0° 

318 
318 
318 
318 
318 
318 
318 
318 
318 
318 
318 
318 
318 
318 
318 
318 


58 
48 
56 
57 
62 
59 
63 
56 
52 
57 
57 
58 
56 
59 
54 
59 


1378 
1373 
1379 
1379 
1383 
1381 
1384 
1379 
1375 
1379 
1379 
1382 
1380 
1380 
1377 
1381 


97 
41 
78 
99 
17 
18 
15 
36 
58 
57 
57 
16 
76 
76 
40 
35 


79324 

27655 
35353 
21071 
20213 
2  0  7  5  9 
24665 
25983 
12968 
11910 
21197 
53414 
58582 
02125 
07746 
25438 


79445 
79204 
79480 
79489 
79627 
79540 
7  9  6  6  9 
79462 
79298 
79471 
79471 
79583 
79522 
79522 
79377 
79548 


869 
884 
047 
747 
729 
740 
823 
850 
576 
554 
749 
426 
535 
349 
466 
838 


T  =309 
s 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

1? 

1  3 

14 

15 

16 


50   0 


32  8 
32  8 
32  8 
32  8 
328 
3  2  8 
32  8 
328 
328 
328 
328 
32  8 
328 
32  8 
328 
328 


21 
2  6 
27 
16 
15 
21 
31 
31 
21 
09 
2  6 
26 
19 
11 
26 
21 


4  5° 

327 
32  7 
327 
327 
327 
327 
32  7 
32  7 
327 
327 
327 
327 
32  7 
327 
32  7 
327 


79 
83 
87 
72 
75 
84 
84 
93 
71 
71 
79 
7  9 
69 
72 
83 
75 


1339 
1342 
1343 
1335 
1336 
1341 
134  1 
1347 
1337 
1333 
1341 
1341 
13  36 
1334 
1342 
1338 


79 
7  4 
98 
93 
61 
36 
36 
48 
14 
65 
36 
5  4 
27 
52 
59 
55 


07511 
19372 
97942 
33687 
00068 
69455 
77899 
80415 
71214 
88105 
44784 
48963 
75752 
92132 
23456 
45635 


2.7 
2  .7 

2  .7 
2  .7 
2  .7 
2  .7 
2  .7 
2  .7 
2  .7 
2.7 


2.7 
2.7 
2  .7 
2.7 


7703 
7833 
7887 
7532 
7562 
7772 
7772 
8042 
7  5  8  5 
7431 
7772 
7780 
7547 
7469 
7826 
7648 


6]] 

866 
406 
173 
450 
793 
974 
029 
980 
193 
412 
502 
077 
279 
954 
430 
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APPENDIX    D 
DATA    REDUCTION    PROGRAM 


Q      -/•  -/\  -X-  -X-  -X-  -X-  -X-  -X-  -X"  -X-  -X-  -X-  "ic  -X-  -5J-  -X-  -X-  -X-  -)c  -X-  -X-  -X-  -X-  *X-  -/?  -X-  -X-  -X-  -X-  -X-  ¥:  -X-  -X-  -)c  -X-  -X-  /c  -X-  -X-  -X-  -X-  -X-  -X-  -X"  -X-  -X-  -X-  -X-  -X-  -X-  -X"  %  "X-  ->c  : 

C         ""-      PROGRAM    TO    CALCULATE    HEAT    TRANSFER    COEFFICIENT    ON    A 
C         -  PINNED    CIRCULAR    DISK 


C 

c  *#*****##*#*#****DEFINITIOn    OF   VARIABLES****************** 

C  A  FACTOR    IN    CHURCHILL    AND    CHU    CORRELATION    USED    HERE 

C  TO    SIMPLIFY    AND    SHORTEN    THE    EXPRESSION 

C  B  LENGTH    OF    THE    SPINE 

C  C  SIMPLIFYING   FACTOR    FOR    CHURCHILL    AND    CHU    CORRELA- 

C  TION 

C  CB  CORRECTION    TO    THE    SPINE    BASE    TEMPERATURE 

C  CO  CORRECTION    TO    THE    SPINE    TIP    TEMPERATURE 

C  D  DIAMETER    OF    SPINE 

C  E  SIMPLIFYING    FACTOR    FOR    THE    CHURCHILL    AND    CHU    COR- 

C  RELATION 

C  G  GRAVITATIONAL    ACCELERATION 

C  GR  GRASHOF    NUMBER 

C  H  HEAT    TRANSFER    COEFFICIENT    (EXPERIMENTAL) 

C  KA  THERMAL    CONDUCTIVITY    OF    THE    FLUID    (AIR) 

C  KP  THERMAL    CONDUCTIVITY    OF    THE    SPINE    MATERIAL 

C  M  =    A    HEAT    TRANSFER    PROPERTY 

C  NU  NUSSELT    NUMBER    (FROM   EXPERIMENTAL    DATA) 

C  NUC  =    NUSSELT    NUMBER    (FROM    CHURCHILL   AND    CHU    CORRELA- 

C  TION) 

C  NUAVG  =    AVERAGE    NUSSELT    NUMBER    (FROM   EXPERIMENTAL    DATA) 

C  NUT  =    TOTAL    OF    NUSSELT    NUMBERS    FROM    16    SPINES    (1RUN) 

C  PHI  =    ROLL    ANGLE    (DEGREES) 

C  PR  PRANDTL   NUMBER    (EXPERIMENTAL) 

C  R  RATIO    OF    TEMPERATURE   EXCESSES    (    =    THB/THO) 

C  RA  RAYLEIGH    NUMBER    (EXPERIMENTAL) 

C  RAT  -    TOTAL    RAYLEIGH   NUMBER    FOR    16    SPINES    (1    RUN) 

C  RAV,  =    AVERAGE    RAYLEIGH    NUMBER    (EXPERIMENTAL) 

C  TB  BASE    TEMPERATURE    (DEGREES    F)    UNCORRECTED 

C  TO  TIP    TEMPERATURE    (DEGREES    F)    UNCORRECTED 

C  TBF  =    CORRECTED    BASE    TEMP.     (DEGREES    F) 

C  TOF  =    CORRECTED    TIP    TEMP.     (DEGREES    F) 

C  TBCD  =    CORRECTED    BASE    TEMPERATURE    (KELVIN) 

C  TO(I)  =    CORRECTED    TIP    TEMPERATURE    (KELVIN) 

C  THB  =    TEMPERATURE    EXCESS    AT    THE    BASE    OF    SPINE 

C  THO  =    TEMPERATURE    EXCESS    AT    THE    TIP    OF    SPINE 
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C     TF      TEMPERATURE  OF  FILM  BETWEEN  SPINE  AND  FLUID 

C     TS      TEMPERATURE  OF  SURROUNDING  FLUID  (AIR) 

C 

c    *x*x#**tt*##*****DECLARATIoN  OF  VARIABLES**************** 

REAL  D , KP , B , TS , R , M , H , PHI , NU , THB , THO , RAT , RAV , NUT , NUAVG , 
REAL  A , C , E , G , TF , VA , ALFA , KA , GR , PR , RA , NUC 

REAL  TB(25),TO(25),CB(25),CO(25) 

INTEGER  I 
C 
C    --------------INITIALIZATION  OF  VARIABLES**************^ 

D  =  0.0127 

B  =  0.0508 

KP  =  120 

G  =  9.807 

TS  =  87. 53 

PHI  =  45 

NUT  =  0 

RAT  =  0 
C 

C    CALCULATE  TS  IN  KELVIN 
C 

TS  =  (TS  -  32)-5/9  +  273.1 5 
C 
C    PRINT  HEADINGS 

WRITE(6,100)  TS,PHI 

WRITE(6,1 50) 
C 

C    READ  IN  AND  CALCULATE  CORRECTED  TEMPERATURES,  DEGREES  F 
C 

DO  50  I  =  1  ,16 

READ(5,200)TB(l) , TO ( I ) , CB ( I ) , CO ( I ) 
TBF  =  TB(I)  +  CBCI) 
TOF  =  TO(I)  +  CO(I) 
C 

C    CONVERT  CORRECTED  BASE  &  TIP  TEMP.  TO  KELVIN 
C 

TB(I)  =  (TBF  -  32)-5/9  +273.15 
TO(I)  =  (TOF  -  32)-5/9  +273.15 
C 

C    CALCULATE  THETA  BASE  AND  THETA  ZERO  (TEMP.  EXCESSES) 
C 

THB  =  TB(I)  -  TS 
THO  =  TO(I)  -  TS 
C 

C    CALCULATE  RATIO  OF  THETA  B  TO  THETA  0  (=cosh(MB)) 
C 

R  =  THB/THO 
C 

C    CALCULATE  M 
C 

M  =  (ALOG(R  +  SQRT( (R--2)-1 ) ) )/B 
C 


C    CALCULATE  HEAT  TRANSFER  COEFFICIENT  H 
C 

H  =  (M**2)*KP*D/4 
C 

C    CALCULATE  FILM  TEMPERATURE 
C 

TF  =  (1  JU.  )*(TB(I)  +  TO(I))  +  (1./2.)*TS 
C 

C    CALCULATE  PROPERTIES  OF  AIR  AT  TF 
C 

KA  =  0.004033  +  (7.30E-5)*TF 

VA  =  (-1.22567  +  (9.350E-3)*TF)*(1 .OE-5) 

ALFA  -  (-1.82333  +  (0. 01 350*TF) )* (1 . OE-5 ) 
C 

C    CALCULATE  NUSSELT  NUMBER  BASED  ON  EXPERIMENTAL  H 
C 

NU  =  H*D/KA 

NUT  =  NUT  +  NU 
C 

C    CALCULATE  GRASHOF,  PRANDTL,  RAYLEIGH,  AND  NUSSELT  NUM- 
C    BERS  BASED  ON  CHURCHILL  &  CHU  CORRELATION 
C 

PR  =  VA/ALFA 

GR  =  (G*(1  .0/TF)*(TF-TS)-;;-D**3)/VA**2 

RA  =  GR  *  PR 

A  =  0.387*RA**(1  .0/6.0) 

C  =  (1.0  +  (0.559/PR)*-;;-  (9.  0/1  6.0  ))**  (8.  0/27.0) 

E  =  A/C 

NUC  =  (0.60  +  E)**2 

RAT  =  RAT  +  RA 
C    PRINT  OUTPUT 
C 

WRITE(6,250)I,TB(I) ,T0 ( I ) ,RA ,NU,NUC ,H 
50    CONTINUE 
C 

RAV  =  RAT/16.0 
NUAVG  =  NUT/1  6.0 
WRITE(6,300)RAV, NUAVG 
C 

STOP 
1 00   FORMAT ('1  ' ,1  5X,  'TS=' ,F6.2,3X, 'PHI=' ,F4.1  ) 
1  50   FORMAT (17X,  'PIN' , 4X,  ' TB ' , 7X,  ' TO ' , 7X,  ' RA '  ,  8X,  ' NU ' , 8X, 

#     'NUC  ,7X, 'H'  ) 
200   FORMATU(F  7.3)  ) 
250   F0RMAT(18X,I2,2X,F6.2,3X,F6.2,2X,F7.2,2X,F8.  5,2X,F8.  5, 

§    3X,F6.3) 
300   FORMAT(1  6X,  'RAV  = ' , F7. 2, 2X, ' NUAVG  =',F6.3) 
C 

END 
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